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Abstract 

In this note we consider a branching Brownian motion (BBM) on R 
in which a particle at spatial position y splits into two at rate f5y 2 , where 
/3 > is a constant. This is a critical breeding rate for BBM in the 
sense that the expected population size blows up in finite time while the 
population size remains finite, almost surely, for all time. We find an 
asymptotic for the almost sure rate of growth of the population. 

AMS 2000 subject classification: 60J80. 
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1 Introduction 

We consider a branching Brownian motion with a quadratic breeding potential. 
Each particle diffuses as a driftless Brownian motion on R, and splits into two 
particles at rate f3y 2 , where {3 > and y is the spatial position of the particle. 
We let N t be the set of particles alive at time t, and then, for each u <E N t , Y u (t) 
is the spatial position of particle u at time t (and for < s < t,Y u (s) is the 
spatial position of the unique ancestor of u alive at time s). We will call this 
process the (f3y 2 \ K)-BBM. 
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It is known that quadratic breeding is a critical rate for population explo- 
sions. If the breeding rate were instead (3\y\ p for p > 2, the population size 
would almost surely explode in a finite time. However, for the (/3y 2 ; R)-BBM 
the expected number of particles blows up in a finite time while the total number 
of particles alive remains finite almost surely, for all time. For p £ [0, 2) the 
expected population size remains finite for all time. See Ito and McKean [3J pp 
200-211] for a proof of these facts using solutions to related differential equa- 
tions. In this note we prove the following result on the almost sure rate of 
growth of \N t \. 

Theorem 1. Suppose that the initial configuration consists of a finite number 
of particles at arbitrary positions in R. Then, almost surely, 



We define Rt := max uS N t Y u (t) to be the right-most particle in the (/3y 2 ; R)- 
BBM. In Harris and Harris [2] it was shown that 



almost surely, and this result is crucial in our proof of Theorem [T] because it 
allows us both some control over the maximum breeding rate in the BBM, and 
also to show that the growth rate of \N t \ is dominated by particles with spatial 
positions near R t . Indeed, a BBM in which every particle has branching rate 
/3i? 2 would see its population grow as in Theorem [TJ and this point will provide 
our upper bound. Furthermore we shall see that, for any S > 0, at all sufficiently 
large times t, a single particle located near R t will single-handedly build during 
time interval [t, t + 5] a progeny so large that it is again of the magnitude given 
by Theorem [1] which will prove the lower bound. 

2 Proof of the growth rate 

To prove Theorem Q] it is sufficient to consider the case of the initial configura- 
tion of particles being a single particle at position iel. 

Upper bound. Since the breeding potential is symmetric about the origin, 
we have from equation (JTJ) that for e > fixed, there exists, almost surely, a 
random time r < oo such that for all u € N t , 



As a consequence of this, the breeding rate of any particle in the population is 
bounded by /3 t := p#WW+*)t for a rj t > t. 

We now introduce a coupled branching process (N t ,t > 0) as follows. For 
each t > it consists in a population of particles {u : u <E N t } whose positions 
are denoted by {Yu{t) : u <E Nt}- Until time r the two processes Nt and Nt 





(1) 



In \Y u (t)\ < (v/2/3 + e)t, for all t > t. 
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coincide and for all t > r we will have N t C N t . Furthermore, we want that 
after time r all particles in N t branch at rate ft. 

In order for this to make sense we must construct (N t ,t > 0) conditionally on 
(N u t > 0). Given (N u t > 0), we let (N t ,t < r) = (N t ,t < r). For J > r, each 
particle u e N t gives birth to an extra particle v G N t \N t at rate ft — (3Y u (t) 2 
(note that if u e iV t then = Y u (t)). Each particle thus created starts an 

independent BBM in (N t ,t > 0) with time-dependent branching rate ft at time 
t. Thus it is clear that all particles in Nt branch at rate ft for t > r. 

It is furthermore clear that 

|JVt| < \N t \ 

for all f > 0. The upshot is that after time r, JV t is a pure birth process, and as 
such is very well studied. 

Consider (Z t ,t > r) a pure birth process starting with a single particle at 
time r with inhomogeneous rate (ft , t > r) , and define 



V t :=inf |s>t: ^ ftdu>t|. 



It is easily seen that (Zv t ,t > 0) is simply a Yule process (i.e. a pure birth 
process where particles split into two at rate 1) and so we know that 

Zy t e~* — > Woo as t — > oo 

almost surely, where Woo has an exponential mean 1 distribution. However, 
since j3 u du — t we see that this implies 

Z t exp ( — / ft ds ) — >• Woo as f — > oo. 



As Nt ioi t > t is composed of the union of the offspring of the |iV r | particles 
present at time r, we can write 

M t := |JVt| exp ( - J ft dsj = ^ Z t W exp ( - y ft dsj 

where the ZM are independent, identically distributed, copies of a pure birth 
process with birth rate ft. By conditioning on the value of |AT T | (which is almost 
surely finite), we see that 

M t -> Moo € (0, oo) as t -> oo 

where, more precisely, Moo is distributed as the sum of |iV T | independent expo- 
nential variables with mean 1. 
Thus t 

In j AT t | — / ftdsWlnMoo as t -> oo 
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and since J /3 S ds — i[exp(ci) — exp(cr)] with c = 2(y / 2/3 + e), it follows that 

In In I Ml , , — 
lim J_*i =2 (v^8 + e), 

t— >oo I 

almost surely. Finally, since e may be arbitrarily small, using the coupling 
constructed above, we have that 

lnln|M| „ 
lim sup — < 2^2/3, 

t— >oo t 

almost surely, as required. 

Lower bound. Fix e > 0. From ([I}, we know that there exists, almost surely, 
a random time r' < oo such that for all t > t' , 

R t > e W-4*. (2) 

Fix 5 > and let t n := nd, n G {0, 1,2,.. .}. We will show that it is sufficient 
to consider only the offspring of the rightmost particles Rt n to prove the lower 
bound for the global population growth rate. There are two steps to the argu- 
ment. First, we will couple the sub-population descended from R tn during the 
time interval [t n , t n+ i] to a BBM in which the expected population size remains 
finite. Second, we will show that, as n tends to infinity, the population sizes of 
the coupled processes grow sufficiently quickly. For this, we use equation ([2} to 
give us a lower bound on the breeding rate in the coupled processes. 

Let Ng 71 ^ be the set of descendants of R tn at time s G [t„, in+i]. As for the 
upper bound we introduce a coupled process (iVj , ,s G [t n ,t n+ i]). This time 
we have that for all s G [t n , t n+1 ],N^ n) C ivj n) C N s . More precisely, (N^ n \s G 
[i n ,i„+i]) is obtained from (Ns , s G [t ni t n+ i]) by cancelling some of the split 
events in Ns 7l \ in the following way: if at time s a particle u G Ns C N^ 
splits in the original process Ns , it also splits in N^ with probability 



2 



G [0,1], 



where 



' /3e 2( v ^-2e)t„ if x > e (VW-2e)t„ 

otherwise. 

(n) 



If the split event is rejected, one of the two offspring in Ns is chosen at random 

— (n) 

to be the one which we keep in iV s . 

The process (Ns n \t n < s < t n+ i) is thus a BBM started at time t n from a 
single particle at position R tn with space-dependent branching rate /3 n . Observe 
that we trivially have N^ C Ns for all s G [t n ,t n+ i\, as announced, since we 
obtain N^ from Ns merely by erasing some particles from N^ n \ 
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We also require another process on the same probability space, coupled to 
(Ns n \s £ [t n ,t n+ i]). This process is denoted (ivj n \ s £ [t n ,t n+ i]), and is 
defined by adding particles to (Ng,s € [i n ,^n+i]) in such a way that every 
particle in breeds at constant rate /j e 2 (VW-^)t n ; irrespective of its spatial 
position. 

More specifically, for s € [t n ,t n+ i], every particle u € gives birth to an 
extra particle v € Ns n \Ns n ^ at rate 

P e L {Y M (s)<e(VW-2e)t n }, 

and each particle thus created initiates an independent BBM in (ivi™' 1 , s € 
[i„, i„+i]) with constant breeding rate f] e 2 (VW-2e)t n , Note that we have Ns ^= 
N { s n) for all s G [t n ,t n+1 ]. 

Now for v € N+ define the event 

AJv) := ( min YJs) < e (VW-2e)t n \ 
I se[t n ,t n+1 ] 1 

and set A n — U e >y(n) A n {v), i.e., the event that there exists a descendant v 

of Rt n (in the modified process N^) such that Yu(s) < e (v /2 ?- 2£ ) t '» for some 
s € [t n ,t n+ i\. Observe that 

{nQ^RQ^CA^ (3) 

We also define the event B n as 

B n :={R tn >e^~^}, 

and let {J-t}t>o be the natural filtration for the (/3y 2 ; R)-BBM. Our aim is to 
show that, almost surely, only finitely many of the events A n occur, from which 
it follows that the populations of the coupled subprocesses are equal to the 
populations of the N^ n ' processes for all sufficiently large n, almost surely. Then 
the final step is to show that the populations 7V( n ) grow sufficiently quickly to 
imply the desired lower bound on the size of the original population. 
We start by writing the event of interest as 

A n = (A n n B n ) U (A n n .B£), 

and we recall from equation ([T]) that P(limsup B c n ) = 0, and hence only finitely 
many of the events A n l~l B^ occur. We now use a standard 'many-to-one' 
argument (see, for example, Hardy and Harris [1 ) to bound the probabilities 
P(A n nB n \T t J. 

Let be the law of a driftless Brownian motion Y started at the point i£l 
(and E x the expectation with respect to this law). Observing that B n e F tn , 
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we have 



l Bn P(A n \T tn ) < l Bn E 



1 B E Rt - 



exp (^%„(F s )d S ); min F s < e ^- 2s ^ 



p (VW~e)t„ 



min Y s < e (^" 2 ^ 
se[o,<5] 



<exp (/3e 2 ^- 2 ^"j) P e 
Using the reflection principle, we obtain that there exists C > such that 

'^-^-(c- 6 *- -1)) 



min Y s <e (V2P-2e)t n \ =2¥ o( Ys 
*e[o,$] / V 



< Cexp 



(v^ff-s)i„^ e - e i„ 



Combining this series of inequalities gives, almost surely, a faster than ex- 
ponentially decaying upper bound on P(A n D -B n |J-( n ), and so we have shown 
that 

Y^P{A n nB n \F tn ) < oo 

n>0 

almost surely. Since A„ n B„ 6 J"t ra+1 , Levy's extension of the Borel-Cantelli 
lemmas (see Williams [U Theorem 12.15]) lets us conclude that, almost surely, 
only finitely many of the events A n n B n occur. Thus only finitely many of the 
events A n occur, almost surely. 

Recalling equation ([3]) we see that, almost surely, there exists a random 
integer n% < oo such that for all n > n%, we have Ns — Ns for all s £ 

[in ) ^n+l] • 

The population size |ivj"' ) | is a Yule process with constant breeding rate 
p e 2(V2P-2e)tn ^ f or s g \t n ,t n+ i\. Hence \N^ J has a geometric distribution 
with parameter 



p n := exp 



,2(v^r-2e)t„ 



If we define 



'In 



exp (e 2 (V^-3e)t„ 



then the probability that the number of particles in N^' is smaller than q n is 

P{\N^l \ < q n ) = 1 - (1 -p n y~ = p n q n + o{p n q n ). 

Using the Borel-Cantelli lemmas again, we see that there exists almost surely a 
random integer < oo such that l-W^J > Qn f° r all n > n 2- 

Finally, for all n > max{ni,ri2} and all t € [t n +i,in+2], we have that there 
are at least as many particles in Nt as there are descendants of Rt„ at time 
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t n +i, which is to say that \N t \ > |A^"|j. Hence, for n sufficiently large, 



-!— ^ > " +1 — > 2v/2^- 7e, 

t t n t 



(4) 



almost surely. (Certainly we must have n > max{ni,ri2}, but we may require 
that n be larger still in order that the multiplicative factor t n /t is close enough 
to 1 for equation (|4|) to hold.) 

We can take e to be arbitrarily small, and so obtain 



almost surely, which completes the proof. 
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